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If QCD is renormalized by minimal subtraction (MS), at higher orders, the strong 
coupling constant a a and the quark masses m q exhibit discontinuities at the flavour 
thresholds, which are controlled by so-called decoupling constants, £ s and f m , 
respectively. Adopting the modified MS (MS) scheme, we derive simple formulae 
which reduce the calculation of £ 9 and £ m to the solution of vacuum integrals. This 
allows us to evaluate f 9 and f m through three loops. We also establish low-energy 
theorems, valid to all orders, which relate the effective couplings of the Higgs boson 
to gluons and light quarks, due to the virtual presence of a heavy quark h, to the 
logarithmic derivatives w.r.t. ra/, of ( g and £ m , respectively. We also consider the 
effective QCD interaction of a CP-odd Higgs boson and verify the Adler-Bardeen 
nonrenormalization theorem at three loops. 



1 Introduction 

It is generally believed that quantum chromodynamics (QCD) is the true the- 
ory of the strong interactions. There are still open questions, concerning the 
origin of confinement or as to why the quark masses m q and the asymptotic 
scale parameter A have the values they happen to have. The answers to these 
questions probably lie outside the scope of perturbative QCD, which forms the 
basis of this presentation. In perturbative QCD, the strong coupling constant 
a s = g 2 /(4ir), where g is the gauge coupling, is small enough to serve as a 
useful expansion parameter, and quarks and gluons may appear as asymptotic 
states of the scattering matrix. 

QCD is a nonabelian Yang-Mills theory based on the gauge group SU(3). 
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In the covariant gauge, the Lagrangian reads 




G^ v — d^G c l — d v G a ^ + gf a G^G1, (1) 

with nf = 6 flavours of quarks ip q (i = 1,2,3), gluons (a = 1,...,8), 
Faddeev-Popov ghosts c a , and gauge parameter £. The generators T a satisfy 
the commutation relations [T a ,T b ] = if abc T c , where f abc are the structure 
constants. 

In the calculation of QCD quantum corrections, one generally encounters, 
among other things, ultraviolet (UV) divergences, which must be regularized 
and removed by renormalization. In quantum electrodynamics, it is natural to 
employ the on-shell renormalization scheme, where the fine-structure constant 
is renormalized in the limit of the photon being on its mass shell. Due to con- 
finement, this limit cannot be taken in QCD, and it is natural to employ the 
most convenient renormalization scheme instead. It has become customary ta 
use dimensional regularizatiom in connection with minimal subtraction (MS)0 
I.e., the integrations over the loop momenta are performed in D = 4 — 2e 
space-time dimensions, introducing a 't Hooft mass scale fj, to keep the (renor- 
malized) coupling constant dimensionless. The poles in e that emerge as UV 
divergences in the physical limit e — > are then combined with the bare (UV- 
divergent) parameters and fields in Eq. (0) so as to render them renormalized 
(UV finite). This is always possible because QCD is a renormalizable theory. 
In the modified MS (MS) scheme^ the specific combination of transcendental 
numbers that always appears along with the poles in e is also subtracted. In the 
following, bare quantities will be denoted by the superscript '0.' Specifically, 
we have 

g° = fJ, e Z g g(iJ,), m° q = Z m m q {^), f-l = Z^(p) - 1), 

1>°f = V^ilftjj), G°>° = ^G>), c°> a = Jz~ 3 c a (v). (2) 

In the MS-like schemes, the renormalization constants Z may be written in 
the simple form 

oo i i 

Z = l + Y.zZ Z ^' (a) 

where a = a s /ir is the renormalized couplant and Zij are numbers. I.e., the 
Z factors do not explicitly depend on dimensionful parameters. In particular, 
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Z m is generic for all q. A crucial advantage of the MS-like schemes is that 
Z g and Z m are £ independent to all orders. This property carries over to a s 
and m q , so that it makes sense to extract these parameters from experimental 
data. 

Z g and Z m carry the full information on how a s and m q run with fi. In 
fact, from the /i independence of g° and rrfi it follows that 



(3(a) 



da 



d\n fi 
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dln/i 2 <iln/! 2 




7m (a) ee = -^Fff = - > : 7«a" +1 - < 4. 



n=0 



The Callan-Symanzik /3 function and the quark mass anomalous dimension 
7 m are universal in the MS-like schemes. Moreover, /3q and /?i are universal in 
the larger class of schemes which have mass- independent functions. In the 
MS-like schemes, the coefficients (3 n and 7„ are known through four loops, i.e., 
n = 3d 

To summarize, the MS-like schemes offer several advantages. They are 
easy to implement in symbolic manipulation programs and are tractable at 
high numbers of loops. Furthermore, a s and m q are £ independent to all 
orders and may thus be regarded as physical observables. The price to pay 
is that heavy quarks do not automatically decouple. However, as will become 
clear in the following, the theoretical ambiguity associated with the matching 
at the flavour thresholds is negligible if higher orders are taken into account. 



2 Decoupling of Heavy Quarks 

The decoupling theorem states that the infrared structure of unbroken, non- 
abelian gauge theories is not. affected by the presence of heavy fields coupled 
to the massless gauge fieldslj As a consequence, a heavy quark h decouples 
from physical observables measured at energy scales fi <C up to terms 
of 0(i-i/rrih). However, the proof of this theorem relies on the use of mass- 
dependent (3 functions. Thus, this theorem does not automatically hold for 
the parameters and fields in MS-like schemes. The standard way out is to 
implement explicit decoupling by using the language of effective field theory. 

As an idealized situation, consider full QCD with ni — rif — \ light quarks q, 
with m q <^ i-i, plus one heavy quark h, with mh 3> /i. The idea is to construct 
an effective theory, QCD', by integrating out the h quark. The parameters 
and fields of the effective theory, which will be denoted by a prime, are related 
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to their counterparts of the full theory by the decoupling relations, 



.g ' = C°<7 , m°' = C° m m° q , - 1 = C 3 °(£° - 1), 

<'* = (C 2 °) 1/2 <\ G2f'° = (C 3 °) 1/2 Gr, c '- = (C 3 °) 1/2 c - a . (5) 

By gauge invariance, the most general form of the effective Lagrangian C 
emerges from Eq. (|l|) by only retaining the light degrees of freedom and reads 

C! (g° s , m°, £°;ipq ,i ,G°^ a , c°>«; C°) = C (gf, m°', £°'; G%", c°'' a ) , (6) 

where £° collectively denotes all decoupling constants of Eq. (||). The latter 
may be derived by imposing the condition that the results for n-particle Green 
functions of light fields in both theories should agree up to terms of 0(n/rrih)- 
As an example, let us consider the g-quark propagator. Up to terms of 
0(fj,/mh), we havaEl 



j[l + TP v (p>)] 



a / <>■<■<»■' WW) = £ w+km * (7) 



where the subscript V reminds us that the self-energy of a massless quark only 
consists of a vector part. Note that Ey(p 2 ) only contains light degrees of free- 
dom, whereas Ey(p 2 ) also receives virtual contributions from the h quark. As 
we arc interested in the limit mh — > oo, we may nullify the external momentum 
p, which entails an enormous technical simplification because then only tad- 
pole integrals have to be considered. In dimensional regularization, one also 
has E^(0) = 0. Thus, we obtain 

C 2 ° = l + £^(0), (8) 

where the superscript h indicates that only diagrams involving closed /i-quark 
loops need to be computed. In a similar fashion, one obtains 

C = i + s|(o) ' C 3° = i + n^(0), c 3 ° = i + n^(o), 

3 = 1 + 1^(0,0), C = ^|^, 0) 

where Eg, Hg, IT, and Tq S c denote the scalar part of the g-quark self-energy, 
the gluon self-energy, the ghost self-energy, and the Gcc vertex function, re- 
spectively. Typical Feynman diagrams contributing to £2; Cm: ds: C 3 : an d Ci 
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Figure 1: Typical Feynman diagrams contributing to (®> Cmi Cjp C3 > an d C' 



are depicted in Fig. [IJ The full set of diagrams .are generated and evaluated 
with the symbolic manipulation packages QGRAIQ and MATAD(3 respectively. 
The renormalized counterparts of Cg and C m , 

C - ^-C° C - —C° (10) 

^g L m 

are found to be UV finite and £ independent and to satisfy the appropriate 
renormalization group equations, which constitutes a strong test. The resulting 
decoupling relations take a particularly simple form if the matching scale is 
chosen to be Hh — mh{^h), namely, 

— = d = 1 + c 2 a 2 + c 3 a 3 , ^ = C™ = 1 + d 2 a 2 + d 3 a 3 , 



in 



'i 



11 564731 82043 >/oS 2633 , 89 

Co = — , c-! = C 3 ni, do — , 

72' 124416 27648 v ' 31104 ' 432' 

2951 In 4 2 ln 2 2 A/nN 407 103,.,, 4 T . (\ 

rf 3 = 1 C(2) C(3) + C(4) Li 4 - 

3 2916 54 9 ^ ' 864 sv ; 72 sv ; 9 V 2 

+ (11664 - 27^)' (U) 

where C and L14 are Riemann's zeta function and the dilogarithm, respec- 
tively. C2 and c?2 were previously calculated^] Three-loop expressions for C2 
and C3j which mayiifi useful for parton model calculations, are available for 
the covariant gaugeE3 

The phenomenological implications of Eqs. (^) and ([ll]) are illustrated 
in Fig. @. For consistency, (n + l)-loop evolution must be accompanied by 
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Figure 2: ^< 5 ) dependence of (a) af ] (M z ) calculated from a { a 4] (M T ) = 0.36 and (b) 

m c (Mz) calculated from [m c = rrS^' ' (ftc) = 1-2 GeV and (M^) = 0.118 with evolution 
at one (dotted), two (dashed), three (dot-dashed), and four (solid) loops and appropriate 
matching. 



n-loop matching. Figure |2|(a) shows how ai 5 \Mz), consistently evaluated 
from ai 4 \M T ) = 0.36 to a given order, depends on the scale measured 
in units of the bottom-quark pole mass Mf, = 4.7 GeV, where the bottom- 
quark threshold is crossed. In Fig. ||(b), the analogous study is performed for 
m [ c\M z ) calculated from (x c = rap (Mc) = 1-2 GeV using a { s 5) {M z ) = 0.118. 
As expected, the dependence on the unphysical scale /i^ is gradually getting 
weaker as we go to higher orders. 



3 Effective Lagrangians and Low-Energy Theorems 

An interesting and perhaps even surprising aspect of £ g and £ TO is that they 
carry the full information about the virtual /i-quark effects on the couplings of 
a CP-even Higgs boson H to gluons and q quarks, respectively. To reveal this 
connection, starting from the bare Yukawa Lagrangian of the full theory, 



•Cyuk — 



V<?=1 



(12) 



where v is the Higgs vacuum expectation value, one integrates out the h quark 
by taking the limit mi —>■ oo and so derives the effective Lagrangian, 

£ v uk = -K E c?o$ = -^^[Oi], (is) 

z=l i=l 
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which is sparine 
dimension fourE 



by a natural basis of composite scalar operators with mass 
The operators, 



o[ = (g; 



,0/\ „/,0/ 



0\ = G°'> a s/^G°'^ v 



0' 5 = (V^c '- b )c '' < \ 



9=1 



(14) 



are only constructed from light degrees of freedom, while all residual depen- 
dence on the h quark resides in the Wilson coefficients Cf. 

The derivation of Cf proceeds similarly to Eq. (^). Considering appropri- 
ate one-particle-irreducible Green functions which contain_a zero-momentum 
insertion of Oh — m h' l Ph' l Ph m ^ nc nml t m h — * 00 > one 

^(-4C? + 2C 4 °) = -±TO(0), C$(C°2 ~ C 3 °) - 1 ~ S°/(0) - ^^"/(O), 

C 2 °c 3 ° = -iflgE^(Q), c 3 °(Q° + c 5 °) = ^1^(0), ffcg = ±a£r^(o,o),(i5) 

with dl = {mfd/dmf), which may be solved for C°. Only 0\ and 2 
contribute to physical observables. They mix under renormalization a£j 



Pi] = 



1 + 2 



Oi-4^1n^)^, [0 / 3 ]=^ 2s (16) 



where the brackets denote the renormalized counterparts. C\ and C 2 are ac- 
cordingly determined from the second equation in Eq. (^). They are dia- 
grammatically calculated through three loopsllllnserting Eqs. (^) and (^) into 
Eqs. ([l5|) , one obtains the low-energy theoremfij 



Ci 



1 91nC 2 



G 2 



1 + 2 



(17) 



which are valid to all orders in a s . Fully exploiting the present knowledge of 
Eq. one may construct the four- loop terms of Q and £ m involving mm 2 
and so obtain C\ and C2 from Eq. ( |l7j ) to one order beyond the diagrammatic 
calculation. The expansions in a = af lf \nh)/'^ reaJ0 



Ci = - — [1 + 2. 7500 a +(9.7951 - 0.6979 ni)a 2 



+ (49.1827- 7.7743 n x - 0.2207 nf)a 3 }, 
C 2 = 1 + 0.2778a 2 + (2.2434 + 0.2454 m)a 3 
+ (2.1800 + 0.3096 m - 0.0100 nf)a 4 . 



(18) 



Having established £y uk , we are able to make higher-order predictions for 
the QCD interactions of a light H boson by just computing massiess diagrams. 
For instance, the H — > gg partial decay width at three loops is found to b<£3 



T(H - gg) = 



G F Mfj 
367iV2 



1 + 17.917a' 



156.808 



5.708 In 



(19) 

where a' = ^^(Mh)/^- The three-loop ©(a^G^ra 2 ) corpptions to T(H — > 
qq), with q — u,d, s, c, b, may also be obtained from Eq. (|l3|)LLj Analogously, the 
QCD interactions of a CP-odd Higgs boson A may be described by an effective 
Lagmiigian involving composite pseudoscalar operators with mass dimension 
fourEil The resulting counterpart of Eq. ( |l9| ) is found to bd 



-> gg) 



16ttV2 



1 + 18.417 a' 



171.544 



r 1 '"t 

5 In — — 4- 

M 2 



(20) 



where a' = (Ma)/7T. As a-by-product of this analysis^ the Adler-Bardeen 
nonrenormalization theorempi which states that the anomaly of the axial- 
vector current is not renormalizcd in QCD, is verified through three loops by 
an explicit diagrammatic calculation. 



4 Comparison with Scale Optimization Procedures 

It is interesting to compare the exact values of the C(a 2 ) corrections in Eqs. ( |l9| ) 
and ( |2C| ) with the estimates one may derive from the knowledge of the 0(a s ) 
correction through the application of well-known scale optimization proce- 
dures, based on Grunberg's concept of fastest apparent convergence (FAC), 
Stevenson's principle of minimal sensitivity (PMS), and the proposal by Brod- 
sky, Lepage, and Mackenzie (BLM) to resum the leading lighl-quark contri- 
bution to the renormalization of the strong coupling constant pi The resulting 
estimates are listed in Table |[ We observe that the sign and the order of 
magnitude is correctly predicted in all cases. 



5 Summary 

A consistent MS description of a s (/j,) and m q (iJ,) with /i evolution through four 
loops and threshold matching through three loops is now available. Effective 
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Table 1: Scale optimization estimates for the 0(a 2 s ) coefficients in Eqs. (hd) and O 





FAC 


PMS 


BLM 


H^gg 


263.3 


263.9 


242.5 


A^gg 


277.6 


278.1 


252.5 



Lagrangians and low-energy theorems are useful tools to treat the hadronic 
decays of light CP-even and CP-odd Higgs bosons through three loops. The 
sign and the order of magnitude of the resulting three-loop corrections are 
correctly predicted by scale optimization procedures. 
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